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Criteria for experimental observation of multi-dimensional optical solitons in media with sat- 
urable refractive nonlinearities are developed. The criteria are applied to actual material param- 
eters (characterizing the cubic self-focusing and quintic self-defocusing nonlinearities, two-photon 
loss, and optical-damage threshold) for various glasses. This way, we identify operation windows for 
soliton formation in these glasses. It is found that two-photon absorption sets stringent limits on the 
windows. We conclude that, while a well-defined window of parameters exists for two-dimensional 
solitons (spatial or spatiotemporal) , for their three-dimensional spatiotemporal counterparts such a 
window does not exist, due to the nonlinear loss in glasses. 

PACS numbers: 42.65.Tg, 42.65.-k 



I. INTRODUCTION 

Solitons are localized wave packets and/or beams that 
result from the balance of the linear and nonlinear re- 
sponses of a physical system. Depending on the physical 
properties of the underlying system, solitons take dif- 
ferent forms. They can be hydrodynamic wave packets, 
such as solitary waves in the ocean 1] and atmosphere 
0. They can also be spin- wave packets, such as mag- 
netic solitons Q . Bose-Einstein condensates provide 
a medium to produce matter- wave solitons 0- Other 
examples of soliton dynamics can be found in a wide 
variety of fields, including astrophysics, plasma physics, 
nuclear physics, and even metabolic biology @, 0, H, , 
among others. Very accurate experiments have been per- 
formed with topological solitons (fluxons) in long Joseph- 
son junctions, including a recent direct observation of 
their macroscopic quantum properties |lOj |. 

Solitons in optics, which are known in their tempo- 
ral, spatial, and spatiotemporal varieties (the latter ones 
being frequently called "light bullets"), constitute, per- 
haps, the most versatile and well-studied (both theoreti- 
cally and experimentally) class of solitons in physics. In 
particular, temporal solitons in optical fibers JdJ have 
recently made a commercial debut in high-speed telecom- 
munications links |TJ |2| ■ It has been pointed out that 
multi-dimensional (multi-D) spatiotemporal optical soli- 
tons can be used in the design of high-speed all-optical 
logic gates and, eventually, in all-optical computation 
and communications systems [T^ . 

The balance of linear and nonlinear dynamical features 
is only the first step in the soliton formation. Securing 
the stability of this balance is the second, equally impor- 
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tant step. A well-known difficulty is that the most com- 
mon optical nonlinearity - the Kerr effect in dielectrics - 
gives rise to soliton solutions which are unstable in more 
than one dimension against the wave collapse, as dis- 
cussed (in particular) in original papers |l4l Il5l | and in 
the review |l6j |. Several mechanisms that can suppress 
the collapse have been investigated. These include satu- 
ration of the Kerr nonlinearity |l7l |. higher-order disper- 
sion or diffraction (also referred to as " non-paraxiality" ) 
|l8l |. multiphoton ionization |l9j| . and stimulated Raman 
scattering (SRS) H^,!^. In particular, importance of 
the multi-photon absorption and SRS for the spatiotem- 
poral self-focusing of light in the Kerr medium was in- 
ferred from experimental data in Ref. p^ . However, 
these mechanisms eventually reduce the intensity and 
cause the pulse to expand in time and space, preclud- 
ing the achievement of multi-dimensional solitons '2^ . 

Different versions of the saturable nonlinearity (which 
implies saturation of the cubic nonlinear susceptibility, 
X^^\ in high-intensity fields) have been studied theo- 
retically in detail. It was shown that both rational 
HI ^ IM |27l and cubic-quintic (CQ) ^3 
variants of the saturation readily support stable two- 
dimensional (2D) and three-dimensional (3D) solitons. A 
difference between them is that the former cannot stabi- 
lize "spinning" solitons with an intrinsic vorticit^ but the 
CQ nonlinearity makes it possible, in the 2D |33 . l33Ll3^ 
and even 3D [33 cases. 

The first observation of a self-trapped beam in a Kerr 
medium was reported by Bjorkholm and Ashkin in 1974 
|36j |. The experiment was done in sodium vapor around 
the D2 transition line, and self-focusing arose from strong 
saturation of the transition (i.e. saturation of the linear 
susceptibility, x^^')- Studies of 2D solitons have made 
rapid progress since the mid-1990's in the study of two 
new nonlinearities featuring saturation. Segev et al. pre- 
dicted that the photorefractive (PR) effect in electro- 
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optic materials could be exploited to create an effective 
saturable nonlinear index of refraction that would sup- 
port solitons [s^l ■ PR, solitons were observed experimen- 
tally soon afterward [s^- In parallel to this, there was 
a resurgence of interest in the so-called cascading non- 
linearity, which is produced by the interaction of two 
or three waves in media with quadratic (x^^^) nonlin- 
ear susceptibility. Both ID and multi-D solitons in the 
quadratic media had been studied theoretically in nu- 
merous works (see reviews [s^ and Stationary 2D 
spatial solitons (in the form of self-supporting cylindri- 
cal beams) were first generated in quadratic media by 
Torruellas et al. Later, Di Trapani et al. observed 
temporal x^^^ solitons 0' ^-^id, finally sp atiotemporal 
solitons were produced by Liu et al. |43 . lil ]. Under 
appropriate conditions, both the PR and cascading non- 
linearities may be modeled as saturable generalizations 
of the Kerr nonlinearity (despite the fact that the PR 
media are, strictly speaking, non-instantaneous, nonlo- 
cal, and anisotropic). However, to date, multi-D solitons 
in true saturable Kerr media have not been observed. 

In this work, we examine the possibility of stabilizing 
solitons (arresting the collapse) in saturable Kerr me- 
dia |l7j , from the perspective of experimental implemen- 
tation. Existing theories provide for parameter regions 
where formation of stable solitons is possible, but neglect 
linear and nonlinear losses, as well as other limitations, 
such as optical damage in high-intensity fields ,45]. First, 
we propose a criterion for acceptable losses, and deter- 
mine the consequences of the loss for the observation of 
soliton-like beams and/or pulses. 

Then, as benchmark saturable Kerr media, we con- 
sider nonlinear glasses. Direct experimental measure- 
ments of the higher-order nonlinearities and nonlinear 
(two-photon) loss in a series of glasses allow us to link the 
theoretical predictions to experimentally relevant values 
of the parameters. As a result, we produce "maps" of the 
experimental-parameter space where 2D and 3D solitons 
can be produced. To our knowledge, this is the first sys- 
tematic analysis of the effects of nonlinear absorption on 
soliton formation in saturable Kerr media. We conclude 
that it should be possible, although challenging, to ex- 
perimentally produce 2D spatial and 2D spatiotemporal 
solitons in homogeneous saturable media. Spatiotempo- 
ral solitons require anomalous group-velocity dispersion 
(GVD). Under conditions relevant to saturation of the 
Kerr nonlinearity, material dispersion is likely to be nor- 
mal. In that case, anomalous GVD might be obtained 
by pulse-tilting e.g. On the other hand, the prospects 
for stabilizing 3D solitons seem poor, even ignoring the 
need for anomalous GVD. This conclusion suggests that 
qualitatively different nonlinearities, such as x^'^\ may 
be more relevant to making light bullets. 

We focus on Gaussian beam profiles, which are proto- 
typical localized solutions. Very recent work has shown 
that nonlinear loss can induce a transition from Gaussian 
to conical waves, which can be stationary and localized 
[4^l47| . The conical waves are very interesting, but rep- 



resent a different regime of wave propagation from that 
considered here. 

The rest of the paper is organized as follows. The theo- 
retical analysis of the necessary conditions for the forma- 
tion of the 2D and 3D solitons is presented in Section 2. 
Results of experimental measurements of the nonlinear 
parameters (cubic and quintic susceptibilities, and two- 
photon loss) in a range of glasses are reported in Section 
3. Final results, in the form of windows in the space 
of physical parameters where the solitons may be gener- 
ated in the experiment, are displayed in Section 4, and 
the paper is concluded by Section 5. 



II. THEORETICAL ANALYSIS OF NECESSARY 
CONDITIONS FOR THE EXISTENCE OF TWO- 
AND THREE-DIMENSIONAL SOLITONS: 
LOSSLESS SYSTEMS 

Evolution of the amplitude E of the electromagnetic 
wave in a lossless Kerr-like medium with anomalous GVD 
obeys the well-known scaled equation pi . l2 ^ l27 l l2 ^ 

iE, + ^{E.,, + Eyy + Eu) + f{I)E^O, (1) 

where z and (x, y) are the propagation and transverse 
coordinates, and t is the reduced temporal variable, and 
/(/) is proportional to the nonlinear correction to the 
refractive index An(/). In the Kerr medium proper, the 
refractive index is n{I) = no + An(/) = no + n2l, which, 
as was mentioned above, gives rise to unstable multi-D 
solitons, including the weakly unstable Townes soliton 
in the 2D case [ig • Upon the propagation, the unstable 
solitons will either spread out or collapse towards a sin- 
gularity, depending on small perturbations added to the 
exact soliton solution. 

Conditions for the soliton formation are usually ex- 
pressed in terms of the normalized energy content, but 
from an experimental point of view it is more rele- 
vant to express the conditions in terms of intensity and 
size (temporal duration and/or transverse width) of the 
pulse/beam. They can also be converted into the dis- 
persion and diffraction lengths, which are characteristics 
of the linear propagation. We transform the results of 
Ref. 15] to estimate the parameters of the 2D and 3D 
solitons in physical units. The transformation is based 
on the fact that the solutions are scalable with the beam 
size. Without losing generality, the estimation also as- 
sumes a Gaussian profile for the solutions. The relations 
between the critical peak intensity necessary for the for- 
mation of the soliton and diffraction length, in SI units 
are: 

I 0-52 {i) (x^) for 2D, 

-^critical ~ \ (2) 

[ 0-79 (g) (x^) for 3D, 
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where Ldiffr = 27rnoWg / Xq is the diffraction length of the 
beam with the waist width wg. Eq. |(2Jl is easy to un- 
derstand for the 2D spatial case. For the 2D spatiotem- 
poral and the 3D case, we have assumed that the light 
bullet experiences anomalous GVD, and has a disper- 
sion length equal to the diffraction length, i.e. we have 
assumed spatiotemporal symmetry for the system, as is 
evident in Eq. Q • Further examination of Eq. shows 
that the beam's power is independent of its size for 2D 
solitons, which is a well-known property of the Townes 
solitons, and the hght-buUet's energy decreases as its size 
decreases in the 3D case 0. 

As said above, two different forms of the saturation of 
the Kerr nonlinearity were previously considered in detail 
theoretically, with An{I) in rational form [21 1211271 l28l| . 

^"<'' = (TT77W ■ 

and CQ (cubic-quintic) [29l|-j33|. 

An{I) = na/ - rul^ , (4) 

with both n2 and positive. Although these two models 
are usually treated separately (and, as mentioned above, 
they produce qualitatively different results for vortex soli- 
tons), they are two approximate forms of the nonlinear 
index for real materials. When the light frequency is 
close to a resonance, Eq. Q describes the system well; 
if the frequency is far away from resonance, Eq. Q is a 
better approximation. When / ^ /gat, Eq. |j2Jl can be 
expanded, becoming equivalent to the CQ model. 



An{I) w n2l - (na/Isat) = n2l ~ nj'^ 



(5) 



with n4 = n2//sat- The two models produce essentially 
different results when the expansion is not valid. 

Critical conditions for the formation of 2D solitons 
in these systems were found numerically by Quiroga- 
Teixeiro et al. (2D), and by Edmundson et al. 
and McLeod et al. [23 for the 3D solitons. From those 
results, we can estimate the necessary experimental pa- 
rameters for both the 2D and 3D case by the transfor- 
mation to physical units. The transformation is based 
on scaling properties of the governing equation JQ). The 
estimate again assumes a Gaussian profile, which yields 



stable 



0.16 {n2/ni) 
1.25(712/774) 



for 2D, 



for 3D. 



(6) 



for the minimum peak intensity needed to launch a stable 
soliton, and 

{0.77Ao-y 770774/712 for 2D, 
(7) 
0.3Ao-y 710714/712 for 3D, 

for the minimum size of the beam. The latter translates 
into the minimum diffraction length, 

3. 68A0714 (710/712)^ for 2D, 
idiffr > { (8) 
0.56Ao7i4 (710/712) for 3D. 



In the derivation of the above equation, we have used 
the result from a CQ model for the 2D case The 
validity of the result can be verified from the fact 
that /stable (714/712) ~ 0.16, which gives an error of 
(/stabie(7i4/7i2))^ ~ 0.025 <C 1 in the expansion of Eq. 
(|5]). This means it is appropriate to use a CQ model to 
determine the boundary where the solitons start to be- 
come stable. On the other hand, the result from a model 
with the form of Eq. Q is used instead for the 3D case 
milil, which yields a result of /stable ~ 1-25/sat and /^at 
can always be expressed in 712 and 714 as described in Eq. 

In general, these results show that the required inten- 
sity decreases with (712/714). This means that a larger 
self-defocusing coefficient 714 makes it easier to arrest col- 
lapse, as expected. On the other hand, a larger 714 also 
makes the beam size larger. This is also understandable, 
since stronger self-defocusing reduces the overall focusing 
effect and makes the beam balanced at a larger size. 



III. THEORETICAL ANALYSIS OF 
NECESSARY CONDITIONS FOR THE 
EXISTENCE OF TWO- AND 
THREE-DIMENSIONAL SOLITONS: THE 
LIMITATIONS DUE TO LOSSES 

Up to this point, the medium was assumed lossless. In 
real materials, saturable nonlinear refraction is accounted 
for by proximity to a certain resonance, which implies in- 
evitable presence of considerable loss. Strictly speaking, 
solitons cannot exist with the loss. Of course, dissipation 
is present in any experiment. The challenge is to build a 
real physical medium which is reasonably close to the the- 
oretical models predicting stable solitons. In particular, 
this implies, as a goal, the identification of materials that 
exhibit the required saturable nonlinear refraction, with 
accompanying losses low enough to allow the observation 
of the essential features of the solitons. Under these con- 
ditions, only soliton- like beams ("quasi-solitons"), rather 
than true solitons, can be produced. Nevertheless, in 
cases where losses are low enough for such quasi-solitons 
to exist (the conditions will be described below) , we refer 
to the objects as "solitons" . 

As candidate optical materials for the soliton genera- 
tion, we focus on glasses, as they offer a number of attrac- 
tive properties jl^jliSIs^]- Their x^^'' susceptibility is, 
generally, well-known, varying from the value of fused sil- 
ica (712 '--^ 3 X 10~-^^cm^/W) up to 1000 times that value. 
The linear and nonlinear susceptibilities of glasses exhibit 
an almost universal behavior that depends largely on the 
reduced photon energy [huj/Eg, where fiuj is the photon 
energy, and Eg is the absorption edge, as defined in Rcfs. 
[isL llg . [53 )■ This results in simple and clear trends 
that can be easily understood. The wide variety of avail- 
able glasses offers fiexibility in the design of experiments. 
Glasses are solid, with uniform isotropic properties that 
make them easy to handle and use. There are recent 
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experimental reports of saturable nonlinearities in some 
chalcogenide glasses [sj • The saturable nonlinearity was 
actually measured with the photon energy above the two- 
photon absorption edge, hence this case is not relevant 
to the pulse propagation, as the loss would be unacccpt- 
ably high. However, these measurements encourage the 
search for situations where the nonlinearity saturation is 
appreciable while the loss is reasonably low. 

It is possible to crudely estimate the conditions that 
will be relevant to soliton formation based on the general 
features of the nonlinearities of glasses. The nonlinear- 
ity of the (2n — l)th order will become significant and 
increase rapidly when the photon energy crosses the n- 
photon resonance. Just as the nonlinear index increases 
rapidly (and is accompanied by two-photon absorption, 
2PA) when hcu/Eg ~ 0.5, we expect n4 to become signif- 
icant (and be accompanied by three-photon absorption, 
SPA) when fiu/Eg ^ O.SS. The requirement that be 
appreciable without excessive 2PA or SPA implies that, 
within the window O.SS < Tilj / Eg < 0.5, the solitons may 
be possible. 

To formulate these conditions in a more accurate form, 
it is necessary to identify a maximum loss level beyond 
which the dynamics deviate significantly from that of a 
lossless system. This issue can be addressed by theoreti- 
cal consideration of quasi-solitons in (weakly) dissipative 
systems. First of all, we fix, as a tolerance limit, an appar- 
ently reasonable value of ^tolerance = 10% peak-intensity 
loss per characteristic (diffraction) length, idiffr- From 
what follows below, it will be clear how altering this def- 
inition may impact the predicted parameter window for 
soliton formation. 

If the loss is produced by 2PA, the corresponding evo- 
lution equation for the peak intensity I{z) is 



dl 
dz 



-/32PA-f^, 



(9) 



where (32PA is the 2PA coefficient. It follows that the loss 
per idiffr (provided that the it is small enough) is A/ w 
— /?2PA-f^-^^diffr- The Substitution of the above definition 
of the tolerance threshold, \AI\/I < -^tolerance, into the 
latter result leads to an upper bound on the intensity: 



I <L 



2PA tolerance 



•^tolerance 
/?2PA^diffr 



(10) 



Notice that the condition implies that the diffraction 
length cannot be too short, hence the upper limit in Eq. 
pO|l cannot be extremely high. 
An analogous result for SPA is 



r < I: 



"^tolerance 



3PA tolerance 



/^SPA-^di 



ffr 



which follows from the evolution equation [cf. Eq. 
^]dl/dz — —P^paI^- However, as will be discussed 
later, in the case relevant to the soliton formation, 2PA 
dominates over SPA. 



On the other hand, within the distance necessary for 
the observation of the soliton, its peak intensity must 
remain above the threshold value to prevent disinte- 
gration of the soliton. Solving Eq. ©, this sets another 
constraint on the intensity: 



1 + Nf32PAloLditfy: 



>L 
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(11) 



where Iq is the initial peak intensity, and TV is the num- 
ber of diffraction lengths required for the experiment. In 
this work, we assume iV = 5, which is sufficient for the 
reliable identification of the soliton ■ Note that 

the condition (|ll|l can never be met if the necessary value 
^stable is too high, 



^stable 



(iV/32PAidiffr)" 



(12) 



In the case of Iq > /max, the overall peak- intensity loss 
with the propagation will be > 50%. We will refer to 
the situation in which /stable > /max as a "loss dominat- 
ing" one, and the opposite as "saturation dominating" , 
since l//stabie and l//max can be viewed, respectively, as 
measures of saturation and loss in the system. When sat- 
uration dominates over the 2PA loss, and hence creation 
of the soliton is possible, Eq. Hll|l can cast into the form 
of a necessary condition for the initial peak power, 



lo > In 



stable 



1 



'stable /^max 



(13) 



The material-damage threshold, /damage, also limits 
the highest possible peak intensity that can be used 
experimentally. Although this threshold depends on 
both the material and pulse duration, we will assume 
/damage — 100 GW/cm^, which is typical for nonhnear 
glasses and pulses with the duration ~ 100 fs. Thus, all 
the above results can be summarized in the form 



/mill < /o < min {/2PA tolerance? /damage} ■ 



(14) 



In a material with known nonlinearity and loss, experi- 
mental observation of the solitons is feasible if the corre- 
sponding window l|14|l exists. 

A somewhat simplified but convenient way to assess 
this is to define a figure of merit (FOM). In the case 

when /damage > /2PA tolerance, 



FOM = log 



/2PA 



tolerance 



log 



^tolerance ( 1 • '^4 2 \ a ~ 



log 



^tolerance 



1.42 



-N) 
for 2D, 

- N 
for 3D. 



(15) 
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If -^'damage IS Smaller than I2PA tolerance, the definition be- 



comes 



FOM EE log 



Id 



log 



amagc 



/damage^ (6.4 - 3.68iV/32PA^ 



2PA is observable even for Ticu/Eg < 0.5 owing to the 
absorption-edge broadening present in all glasses. 

Typical experimental traces obtained from AS2S3 are 
shown in the insets of Fig. ^ along with the theoretical 
fits. The intensity dependence of SRTBC signal magni- 
tude and normalized nonlinear absorption signal magni- 
tude are shown in Fig. ^ The dotted curves in both 



log 



^damagcif (o.8 - 0MN(3; 



2PA- 



for 2D, 



Anrin 



(16) 



for 3D. 



The FOM is a measure of the range between the mini- 
mum required and maximum allowed values of the peak 
intensity. Of course, it must be positive, and the larger 
the FOM, the better the chance to observe solitons. 

It seems to be commonly accepted that a larger quintic 
self-defocusing coefficient 77.4 is always desirable, but the 
above results show that this is not always true. From the 
FOM we can see that a larger 714 is better in the sense 
that it reduces the lower threshold /min, helping to secure 
the positiveness of the FOM (|16|1 . However, as soon as 
^min is low cnough that the damage threshold no longer 
poses a problem, Eq. (|15|l shows that larger does not 
help, and the loss factor /32PA dominates. One can under- 
stand this, noticing that, although larger 714 reduces /min, 
at the same time it increases the beam's width and makes 
the needed experimental propagation length longer, as is 
clearly shown by Eq. (jSJ. In turn, more loss accumu- 
lates due to a longer propagation length, which offsets 
the benefit of a lower /min. 



IV. MEASUREMENTS OF NONLINEAR 
PARAMETERS OF GLASSES 

The eventual objective is to answer the following ques- 
tion: for a given category of materials (such as glasses), 
with known nonlinear, loss, and damage characteristics, 
does there exist a combination of material and wave- 
length such that solitons can be observed? To this end, 
we have measured the nonlinearity in a series of glasses 
with 100-fs pulses from a Ti:sapphire regenerative ampli- 
fier with center wavelength at 790 nm. Sapphire is used 
(it has hui/Eg ^ 0.25 in this case) as a reference material 
with minimum nonlinearity. Although fused silica can 
also be used for this purpose, sapphire's higher damage 
threshold allows us to measure at higher intensities. 

We measured several glasses, including: SF59 (with 
huj/Eg ~ 0.5), La-Ga-S(with Hco/Eg ~ 0.56), and AS2S3 
(with hui/Eg ~ 0.75). To determine the effective x^'^-' and 
X^^) susceptibilities, spectrally resolved two-beam cou- 
pling (SRTBC) was used We extended the appli- 
cation of this method to take into account both higher- 
order nonlinearities and strong signals |53j |. In general. 
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FIG. 1: Intensity dependence of (a) SRTBC signal magnitude 
(normalized peak-valley transmission difference) and (b) non- 
linear absorption signal of AS2S3. The saturation of the cubic 
nonlinearity is evident. Higher-order nonlinearities, such as 
X^^\ can be estimated from the deviation from x^^'^ ■ Insets 
show examples of SRTBC and nonlinear absorption traces 
(symbols) along with the best-fit theoretical curves (solid 
lines). The time delay (on the horizontal axis) is given in 
units of the pulse duration (FWHM). 

panels are predictions for the pure x'''^'' nonlinearity. The 
deviation of the experimental points from these curves 
evidences the saturation of the nonlinearity. Postulat- 
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ing the presence of the x^^^ self-defocusing nonhnearity 
provides for good agreement with the experiments. Sim- 
ilar results were produced by all four samples used in the 
measurements; in particular, in all the cases the sign of 
the real part of x^^-* turns out to be opposite to that of 
X^^\ i-e., the quintic nonlinearity is self-defocusing in- 
deed. The measured x^^-* coefficients are consistent with 
previously reported values 48, 51, 54]. 

From these results, we also observe that higher-order 
nonlinearities become more important as the optical fre- 
quency approaches a resonance, as expected on physical 
grounds. The x^^^ part of the nonlinearity is most signifi- 
cant for AS2S3, while for sapphire it is below the detection 
threshold. 



V. STABILITY WINDOWS FOR THE TWO- 
AND THREE-DIMENSIONAL SOLITONS 

The measurements provide the information needed to 
construct the window for the soliton formation. The re- 
sults for 2D case are shown graphically in Fig. |2| The in- 
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FIG. 2: The operation window for the 2D solitons, as pre- 
dicted on the basis of the experimentally-measured character- 
istics of the glass. The hatched area is the window neglecting 
the loss. The shaded area is the dramatically reduced (but 
definitely existing) window found with the loss taken into ac- 
count. 



To demonstrate the dramatic effect of the loss, we also 
plot the window for the (unrealistic) case when loss is 
completely neglected (the hatched area). In the absence 
of loss, the window is very large and the FOM increases 
monotonically with the reduced photon energy. The 
shaded area is the window remaining after inclusion of 
the loss. It is greatly reduced compared to the lossless 
case, and the best FOM is obtained near hui/Eg ~ 0.35. 
From this diagram, we conclude that, while the satura- 
tion of the nonlinearity is definitely necessary to stabilize 
the soliton, major restrictions on the window are imposed 
by the loss. 

From the above rough estimation that were based on 
the band-edge arguments, one might expect that 3PA 
would further curtail the window, when the 2PA effects 
are weak (which is the case exactly inside the predicted 
window). However, n2 and 2PA have been observed in 
glasses for the reduced photon energy as low as ~ 0.35 
|49| . due to the fact that the band edge in glasses ex- 
tends well below the nominal value. Since significant 
2PA remains in this region, 3PA may be neglected in- 
deed. Hence, 2PA presents the fundamental limitation 
to observing solitons in these media [as quantified by the 
FOM in Eqs. (HHJl and ifTB)) ]. 

The results of the analysis for the 3D solitons are sum- 
marized in Fig. 121 Note that another major issue in this 
case is the requirement of anomalous GVD. This require- 
ment is neglected here (addition of it will only further 
constrain the window, which does not really exist even 
without that, see below). From Fig. |31 we observe that. 
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tensity limitations are plotted on the diagram against the 
reduced photon energy. The parameter space can be di- 
vided into two regions which were defined above, viz., the 
saturation-dominating and absorption-dominating ones, 
with the boundary between then determined by Eq. (|12|l . 



FIG. 3: The operation window for 3D solitons. The meaning 
of the hatched area is the same as in the 2D case, i.e., it 
shows the window obtained neglecting the loss. When the 
loss is taken into account, the window vanishes completely. 
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even in the lossless case, the window (hatched area) is sig- 
nificantly smaller than in the 2D case. This is expected, 
because collapse is stronger in 3D than 2D 16]. As in 
the 2D case, the loss again is a major concern for per- 
forming experiments. The most important inference is 
that the window closes up completely when loss is taken 
into account. Thus, it appears that loss will preclude the 
creation of 3D solitons in glasses, while leaving room for 
the 2D solitons. 

Our overall conclusion is that a challenge is to perform 
experimental studies of 2D solitons in saturable Kerr me- 
dia. Both spatial and spatiotemporal solitons are possi- 
ble to be produced experimentally. Among these two, 
the 2D spatiotemporal case is more complicated since it 
requires anomalous GVD. In general, this will naturally 
constrain the window further. On the other hand, in 
this case tilted-pulse techniques could be used to obtain 
anomalous GVD. It is also possible to use a planar waveg- 
uide to perform 2D spatiotemporal soliton experiments. 

Of course, the predicted window depends on the 
assumed parameters (such as the damage threshold) 
and criteria (such as the 10% loss per diffraction 
length) .Variations in these parameters will naturally im- 
pact the window, and our analysis provides the guide- 
lines for searching for the most favorable materials and 
wavelength. A next natural step is to perform numerical 
simulations of the pulse propagation with the parameters 
selected in the present work. It is conceivable that the 
window for 3D solitons would finally open through vari- 
ations of material parameters. In that case, one would 
still have to find an overlap of the resulting window with 
the condition that the GVD must be anomalous. More 
generally, non-glass materials may be tried to improve 
the possibilities for the experiment. 

VI. CONCLUSION 

We have developed criteria for experimental obser- 
vation of multi-dimensional solitons - spatial and spa- 



tiotemporal 2D solitons and spatiotemporal 3D ones. 
Using these criteria and measured properties of nonlin- 
ear glasses within a range of reduced photon energies, 
we have shown that the loss that accompanies higher- 
order nonlinearities (which are tantamount to saturation 
of the cubic nonlinearity) will set very stringent limits on 
the material parameters appropriate for the experiment. 
While loss was thus far neglected in theoretical treat- 
ments of multi-dimensional solitons, this work motivates 
more systematic studies of the soliton-like propagation in 
lossy media. 

The criteria developed in this paper can also be ap- 
plied, as an assessment tool, to materials other that 
glasses. More generally, the same rationale used for ob- 
taining the relevant boundaries in this paper can also be 
used in systems other than optical ones. In these cases 
the specific mathematical forms of the boundaries will be 
different. In any case, the analysis presented here sug- 
gests that there is a small but apparently usable window 
of parameters in which 2D solitons can be generated, and 
work is underway to address this possibility. On the other 
hand, the prospects for generating 3D solitons in glasses 
are quite poor. 
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